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Final  Report  for  AFOSR  88-0015 
for  the  Period 

October  15,1988  to  October  15,1989 
Modeling  and  Control  of  Flexible  Vehicles  in  Space 


Over  the  period  15  September  1987  to  15  October  1989  several  aspects  of  the 
modeling  and  control  of  flexible  space  structures  have  been  investigated: 

1.  The  reduction  of  modal  data  to  forms  useful  for  control  system  designers. 

2.  Design  of  shape  control  systems  for  two  and  three  dimensional  structures 
with  many  vibration  modes  by  decoupling  into  smaller  subsystems,  each 
with  fewer  controls  and  outputs. 


Transformation  of  Modal  Data  from  Finite  Element  Codes  to  a  Format  for 

Control  Design 

The  efficient  reduction  of  modal  data  obtained  from  finite  element  codes  to  a 
form  useful  for  control  system  designers  was  realized  through  the  use  of  a 
program  we  developed  called  TRANSGEN.  This  program  is  useful  for 
producing  state  space  models  of  two  or  three  dimensional  structures  having 
many  sensors  and  actuators.TRANSGEN  accepts  as  input  the  output  of  PAL,  a 
personal  computer  version  of  NASTRAN.  Codes  of  this  type  exist  for 
manipulating  NASTRAN  output  [1],  but  we  believe  this  is  the  first  such  code 
produced  for  PAL  output. 

While  the  generation  of  state  space  models  for  flexible  beams  or  other  one 
dimensional  structures  can  be  performed  by  using  standard  techniques,  such 
as  the  assumed  modes  method,  more  complex  multi-dimensional  structures 
must  be  analyzed  using  finite  element  methods.  Typically  the  designer  of  a 
control  system  for  such  a  structure  is  given  a  modal  form  of  the  system 
equations  of  motion  and  must  determine: 

(1)  Which  sensor  and  actuator  pairs  should  be  chosen  for  control  design. 
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(2)  Which  of  the  many  modes  to  keep  in  the  system  description 

(3)  The  relative  importance  of  each  mode  for  each  individual  sensor/actuator 
transfer  function. 

By  obtaining  the  pole/residue,  pole/zero,  and  state  space  forms  of  the  system 
description  automatically,  TRANSGEN  allows  the  designer  to  experiment 
with  different  linear  or  rotational  sensor  and  actuator  placements  and 
different  model  order  reductions.  A  complete  description  of  TRANSGEN, 
including  operating  instructions  and  examples,  is  provided  in  Appendix  1. 

Two-Dimensional  Flexible  Structure  Control  Design  by  Input/Output 
Decoupling  into  Smaller  Subsystems 


Modelling 

One  of  the  salient  features  of  many  flexible  space  structures  is  a  considerable 
number  of  modes  with  closely  spaced  (or  even  repeated)  eigenvalues  with  low 
damping.  This  arises  from  symmetries  or  periodicity  in  the  geometric  shape 
of  the  structure. 

As  a  representative  space  structure  with  these  features,  we  chose  a  flat, 
circular  plate  whose  control  had  been  investigated  experimentally  by  previous 
researchers  [2].  By  modeling  the  plate  with  lumped  masses  representing 
actuator/sensor  packages  at  symmetric  locations  about  the  x  and  y  axes,  we 
were  able  to  retain  symmetry  in  the  mode  shapes  while  increasing  the 
realism  of  the  model.  The  masses  caused  some  small  symmetric  warping  of 
the  mode  shapes  compared  to  the  mode  shapes  of  the  plate  alone,  as  described 
below. 


The  finite  element  discretization  used  in  PAL  ,  as  well  as  the  placement  of  the 
sensor/actuator  lumped  masses,  is  depicted  in  Fig  1.  Figure  2  shows  the  mode 
shapes  for  the  circular  plate  if  no  lumped  masses  are  added,  while  Figure  3 
shows  the  mode  shapes  for  the  plate  with  lumped  masses.As  seen  in  Figure  2, 
nodal  (or  zero  displacement  )  lines  lie  at  various  increments  of  30  and  45 
degrees  around  the  plate.  A  translation  sensor  or  force  actuator  placed 
anywhere  along  these  lines  will  cause  that  mode  to  be  either  unobservable  or 
uncontrollable,  respectively,  from  that  sensor  or  actuator.  Hence  we  decided  to 
place  4  sensor/actuator  pairs  symmetrically  about  the  x  and  y  axes  as  seen  in 
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Figure  1.  Each  pair  was  modelled  as  a  lumped  mass  of  l/100th  of  the  plate 
mass.  An  additional  lumped  mass  was  included  at  the  center  of  the  plate  to 
approximate  the  mass  of  a  rigid  feed  tower. 

From  Figure  3  we  see  that  modes  number  4  &  8  show  considerable  symmetric 
warping,  but  translation  sensors  and  actuators  at  15°  will  still  be  able  to 
observe  and  control  all  the  modes  depicted  since  none  have  nodal  lines  at  those 
locations.  Placing  the  sensor/actuator  pairs  at  the  edge  of  the  plate  provides 
maximum  observability/controllabillity  since  most  mode  shapes  have 
maximum  deflection  at  the  edge. 

TRANSGEN  was  used  to  create  a  16th  order  state  space  model  of  the  plate 
assuming  4  sensor/actuator  pairs  placed  at  ±  15  degrees  from  the  positive  and 
negative  y  axis.  The  resulting  system  F,G,  and  H  matrices  are  shown  in 
Table  1. 

Control  Design 

By  examining  the  mode  shapes  it  was  seen  that  they  were  separable  into  4 
groups:  symmetric  about  the  x  axis,  symmetric  about  the  y  axis,  symmetric 
about  the  x  and  y  axes,  and  anti-symmetric  about  the  x  and  y  axes.  The 
groupings  are  shown  in  Figure  4  .  By  forming  linear  combininations  of  the 
inputs  and  outputs  of  the  system,  the  G  and  H  matrices  were  decoupled  into 
four  subsystems,  as  shown  in  Figure  5.This  proceedure  is  similar  to  one 
described  in  [3]. 

The  symmetry  to  subsystem  correspondence  is  shown  in  Figure  4.  Figure  4 
also  indicates  whether  -options  about,  or  translations  of,  the  x  and  y  axes 
occur  for  each  subsystem  ,T'his  information  is  useful  since  a  quarter  plate 
could  be  modeled  with  thes^.  sets  of  boundary  conditions  in  order  to  produce  the 
same  modal  displacement  and  frequency  information.  Doing  so  greatly 
reduces  the  size  of  the  finite  element  model  and  thus  allows  faster  analysis. 

Each  subsystem  is  SISO  (single  input  single  output)  and  has  the  properties  of  a 
collocated  system  ,  i.e.  alternating  poles  and  zeros  in  the  transfer  function. 
Simple  lead-controllers  could  then  be  used  to  damp  each  mode.r  Root  loci 
versus  lead  gain  are  illustrated  in  Figures  6  through  9  for  the  four  subsystems. 
The  unmodeled,  higher  frequency  modes  of  each  subsystem  could  also  be 
damped  by  the  lead  controllers  since  they  would  follow  an  alternating  pole- 
zero  pattern.  The  resulting  compensator  for  the  entire  plate  was  therefore 
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fourth  order. 


The  performance  of  the  closed  loop  system  was  evaluated  by  simulating 
impulse  responses  for  a  variety  of  disturbances.  Figure  10a  shows  the  rotation 
about  the  x-axis  at  the  center  of  the  plate  (not  sensed,  but  constructed  from  the 
mode  responses)  for  an  impulsive  disturbance  torque  about  the  x-axis.  Figure 
10b  shows  the  output  of  sensor  1  (z  translation),  and  figures  10c  and  d  show  the 
control  force  from  actuators  1  and  2  (also  in  the  z  direction).  Sensors  2,  3,  and  4 
showed  responses  either  equal  to  the  response  of  sensor  1  or  just  the  inverse  of 
it.  Actuators  3  and  4  responded  exactly  opposite  to  actuators  1  and  2. 

The  plots  in  Figure  ll(a-d)  are  analogous  to  those  in  Figure  10  but  for  a  y-axis 
impulsive  disturbance  torque.  Figures  12  show  the  results  of  a  force  impulse  at 
the  location  of  actuator  1. 

Impulse  responses  for  the  total  closed  loop  system  showed  well-damped 
behavior  and  it  didn't  appear  that  any  of  the  controllers  were  "fighting"  each 
other,  as  might  occur  with  poorly  designed  controllers  for  multi-input  multi¬ 
output  systems. 


Conclusion 

Having  gained  better  understanding  of  the  methods  involved  in  modeling 
large,  two-dimensional  space  structures  via  the  finite  element  approach,  and 
some  of  the  simplifications  inherent  in  the  control  system  design  by  utilizing 
system  symmetry,  we  believed  a  logical  extension  of  our  work  would  be  to 
incorporate  simutaneous  structural  and  control  system  design.  This  is 
thought  to  be  quite  useful  since  a  structure  can  be  designed  to  be  much  lighter 
if  an  active  control  system  is  present  to  minimize  structural  deflections,  and 
therefore  strains,  in  response  to  dynamic  loads  from  disturbances,  slew 
maneuvers,  etc.  Our  research  was  proceeding  in  this  direction  when  AFOSR 
funding  was  stopped. 
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Introduction 

Given  the  modal  output  of  an  undamped  vibratory  system  from  MSC  PAL, 
TRANSGEN  will  calculate  transfer  functions  for  multi-input  multi-output 
(MIMO)  systems,  the  corresponding  pole/residue  form,  and  the  state  space 
form  of  the  system  equations  for  any  selected  input  and  output  locations. 
The  required  PAL  data  are  the  modal  displacements  from  a  normal  modes 
analysis  and  the  modal  masses  from  a  transient  response  analysis 
(unfortunately,  the  PAL  normal  modes  analysis  does  not  have  an  option  to 
print  out  the  modal  masses).  The  required  user  input  is  simply  the  number 
of  nodes  in  the  PAL  model,  the  number  of  modes  to  include  in  calculating 
the  desired  outputs,  the  number  of  actuators  and  sensors,  and  their 
locations  and  orientations.  The  user  also  has  the  option  of  printing  out  any 
combination  of  output  forms  (i.e.  pole-residue,  pole-zero,  and/or  state  space 
models). 

Listed  below  are  instructions  for  using  TRANSGEN,  the  relevant  theory 
behind  the  algorithms  used  in  TRANSGEN,  an  example  problem,  and  the 
FORTRAN  program  listing. 


Producing  PAL  Modal  Data 

After  building  the  PAL  model  file,  both  normal  modes  analysis  and 
transient  analysis  (under  the  PAL  dynamics  menu)  must  be  run  in  order  to 
produce  the  modal  displacements  and  mass  files. 

When  running  the  normal  modes  analysis, the  user  must  specify 
"jacobi"  in  the  input  file  if  the  structure  has  any  possible  rigid  body 
motions.  The  user  may  also  specify  how  many  rigid  body  modes  will  exist  in 
this  file  by  including  the  line  "RIGID-BODY  MODES  n"  where  n  is  the 
number  of  rigid  body  modes  expected.  The  addition  of  this  line  causes  the 
frequencies  of  the  first  n  modes  to  be  exactly  zero  (instead  of  a  very  small 
number)  and  has  little  effect  on  the  mode  shapes.  When  prompted  by  PAL 
the  user  should  choose  the  "print  modal  displacements"  option  from  the 
output  choices,  and  enter  an  appropriate  file  name.  PAL  will  also  ask  for 
the  number  of  modes  to  be  calculated.  Since  dynamic  analysis  can  have  a 
maximum  of  100  active  degrees  of  freedom,  it  is  best  to  choose  less  than  30 
modes  since  mode  shape  accuracy  degrades  above  that. 

As  mentioned  above,  the  modal  masses  are  also  required  for  transfer 
function  calculation  but  cannot  be  output  from  the  normal  modes  analysis 
section  of  PAL.  However,  they  can  be  obtained  from  transient  analysis  by 
using  the  following  input  file, 
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TITLE  EXAMPLE 
JACOBI 

USE  FIRST  n  MODES 
TIME  RANGE  1,1,1 
PRINT  MODAL  EQUATIONS 
SOLVE 

Here,  n  signifies  how  many  modal  masses  are  required  to  be  output,  and 
jacobi  has  been  selected  to  signify  that  rigid  body  modes  will  occur.  It  turns 
out  that  the  number  of  modes  used  has  no  effect  on  the  computation  time. 
When  the  user  is  prompted  for  desired  outputs  he  should  not  select  any  as 
they  would  then  overwrite  the  modal  equations.  The  modal  equations  are 
really  just  the  elements  of  the  diagonal  modal  mass  and  stiffness  matrices. 
An  output  file  name  should  still  be  entered.  As  the  bulk  of  the  time  spent  in 
performing  a  transient  analysis  is  for  recalculating  the  eigenvalues  and 
eigenvectors  (frequencies  and  mode  shapes),  this  run  will  take  as  long  as 
the  normal  modes  analysis,  but  there  is  no  alternative  . 

Running  TRANSGEN 

Once  the  modal  displacement  and  mass  files  are  obtained  TRANSGEN 
can  be  run.  The  user  has  the  option  of  inputting  all  information  by  hand  or 
having  it  read  from  an  input  file,  whose  structure  is  shown  below. 

freebeam.mode 

freebeam.mass 

8.6 

2,2 

1 .3 

6.3 
1  ,3 
6,3 
1.1.1 

The  first  line  contains  the  name  of  the  modal  displacements  output  file.  The 
second  line  has  the  name  of  the  modal  mass  file.  On  the  third  line  is  the 
number  of  nodes  in  the  structure  and  the  number  of  modes  desired  to  be 
used  in  the  subsequent  calculations.  The  number  of  modes  must  obviously 
be  less  than  or  equal  to  the  number  of  modes  for  which  the  user  has 
obtained  modal  displacments  and  masses.  On  the  next  line  is  the  number  of 
actuators  and  the  number  of  sensors,  each  of  which  must  be  <  6  (this  limit 
can  be  increased  by  the  user  by  altering  the  parameter  "ma"  in  the  fortran 
source  file).  The  following  lines  contain  information  about  each  actuator 
and  sensor  starting  with  actuator  number  1  and  ending  with  the  last 
sensor.  Each  line  contains  two  numbers;  the  node  at  which  the  sensor  or 
actuator  is  located,  and  a  number  between  1  and  6  specifying  the  degree  of 
freedom  in  which  it  is  "sensitive".  The  number  is  chosen  according  to  the 
following  scheme, 
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1= 

x  translation 

2= 

y  translation 

3= 

z  translation 

4= 

x  rotation 

5= 

y  rotation 

6= 

z  rotation. 

The  last  line  contains  three  switches  where  l=on,  0=off,  for  each  of  the 
three  output  options,  namely,  residues  output,  transfer  function  zeros 
output,  and  state  space  form  output. 

TRANSGEN  Output 

The  first  section  of  the  TRANSGEN  output  contains  the  zeros  of  the 
transfer  function  numerator  polynomial  matrix  (see  the  section  on  theory 
below  for  more  details  about  this).  This  matrix  is  size  noxni  where  no= 
number  of  sensors  and  ni=  number  of  actuators.  Its  form  is  as  seen  below, 
where  nn(s)  is  the  numerator  of  the  transfer  function  from  actuator  1  to 

sensor  1,  ni;^(s)  is  the  numerator  of  the  transfer  function  from  actuator  2  to 
sensor  1,  and  so  on.  All  transfer  functions  have  the  same  denominator  d(s). 


'y  /s)  1 

n ,  1  (S)  n,2<S> 

'^(s)  ' 

y2(s) 

= 

n21(S)  n22<s>- 

u2(s) 

u  (s) 

.  nr  'j 

d(s) 


d(s)  =  n(s2  +  co.2) 

i 

co^pole  frequencies  =  normal  mode  frequencies 

The  real  and  imaginary  parts  of  each  root  are  printed  since  if  sensor 
and  actuator  are  separated  (non-collocated)  then  zeros  may  exist  which  are 
not  on  the  imaginary  axis  in  the  s-plane  (s  is  the  Laplace  transform 
variable). 

In  all  collocated  cases  all  the  zeros  will  lie  on  the  imaginary  axis  and 
thus  should  have  real  part  equal  to  zero.  Due  to  numerical  inaccuracies 
and  truncation  effects  TRANSGEN  will  yield  zeros  with  non-zero  real 
parts,  but  they  will  be  small  compared  to  the  imaginary  part  and  can  be 
ignored. 

For  non-collocated  cases  the  user  must  exercise  his  own  judgement  for 
neglecting  the  real  part  of  certain  zeros  since  they  may  sometimes  lie  on  the 
real  axis.  If  a  zero  is  found  which  has  an  imaginary  part  within  1  percent 
of  a  pole  (normal  mode)  frequency  the  message  "possible  pole  zero 
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cancellation"  will  appear  following  the  zero's  imaginary  part.  The  user  can 
then  examine  the  pole  frequencies  to  determine  if  he  wishes  to  actually 
consider  them  cancelled. 

The  next  section  contains  a  list  of  the  pole  frequencies,  or  roots  of  d(s) 
(these  are  just  the  normal  mode  frequencies). 

The  third  section  contains  residue  matrices  for  the  first  n  modes  where 
n  is  the  number  of  modes  the  user  wanted  included  in  the  calculations. 
Each  residue  matrix  is  size  no  x  ni. 

The  last  section  of  the  output  contains  the  state  space  model  G  and  H 
matrices  for  the  state  space  equations, 

n+o)D2n  =  Gmu 

y  =  Hmn 

with  q=  vector  of  modal  deflection  amplitudes.  The  dimension  of  Gm  is  n  x 

ni,  and  the  dimension  of  Hm  is  no  x  n.  The  coD2  matrix  is  the  diagonal 

matrix  of  modal  frequencies  squared.  The  user  should  note  that  these  are 
second  order  equations  which  can  be  transformed  into  first  order  equations 
by  the  usual  methods. 
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Theory 

In  this  section  the  algorithms  for  calculating  the  MIMO  residues, 
transfer  functions  and  pole-zero  form  from  the  system  eigenvalues  and 
eigenvectors  (mode  frequencies  and  shapes)  are  explained.  These 
algorithms  are  explained  in  more  depth  by  Chu  and  Plescia  [A1J. 

The  "residues’1  of  a  residue-pole  form  of  a  transfer  function  are  the 
numerators  of  a  partial  fraction  expansion  of  the  pole-zero  form  .  However, 
the  residues  can  be  calculated  without  first  forming  the  pole-zero  form  as 
follows; 

The  system  equations  from  a  finite  element  analysis  can  be  expressed  as 
a  set  of  coupled  second-order  equations: 

M  x  +  K  x  =  Gu  . 

x  =  6n  x  1  nodal  displacements  vector, 
n  =  number  of  nodes  in  structure, 

M  a  mass  matrix  (square), 

K  =  stiffness  matrix  (square), 
u  =ni  x  1  control  input  vector, 

G  =  input  distribution  matrix, 

Then,  using  the  PAL  generated  eigenvector  matrix,  P,  the  system  is 
diagonalized  to, 


MDf|+KDn=PTGu  , 
using  the  transformation, 

x=Pt]  , 

with 

md=ptmp, 

Kd=PTKP,  Mp  KQ diagonal  , 

Md  Kd  are  modal  mass  and  stifness  matrices 

formed  in  PAL  "transient  response  analysis", 
t|  =  m  x  1  modal  amplitudes  vector  , 
m  a  number  of  modes  user  selects  to  be  used  . 


Then  with 


y  =  Hx  =  HPi] , 

it  is  simple  to  form  the  transfer  function  matrix, 
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y(s)  =  h 


1 

i  =1 


\ 

Gu(s)  , 

) 


where  Pi  is  the  i'th  column  of  P  (i  th  mode  shape),  mj  is  the  ith  modal 
mass(the  ith  diagonal  member  of  the  Mp  matrix)  ,  and  Wj  is  the  i  th  modal 
frequency  =sqrt(ki/mi).  Since  H  and  G  are  constants  they  can  be  taken 
inside  the  summation  to  greatly  reduce  the  size  of  the  PjP^  matrix  since 
there  will,  in  general,  be  fewer  inputs  and  outputs  then  the  number  of 
modes  (m).  The  term  HPiPiTG/mi  is  then  the  residue  matrix  for  the  i  th 


mode  (hereafter  called  Rp  whose  frequency  is  coj.  Thus 


y(s)  = 


u(s)  . 


Although  P  is  a  large  matrix  (size  6n  x  m),  the  computation  of  Rj  is  not  as 

difficult  as  it  seems  since  the  number  of  inputs(actuators)  and  outputs 
(sensors)  is  much  smaller  than  the  number  of  nodes  in  the  structure  so  the 
G  and  H  matrices  will  be  sparse  (i.e.  have  lots  of  zeros).  This  greatly 

reduces  the  task  of  multiplying  HPj  (and  P^G)  since  we  have  only  to  select 

the  components  of  eigenvector  i  corresponding  to  the  degrees  of  freedom  in 
which  the  sensors  are  sensitive  (or  the  actuators  provide  force  or  torque). 
Note,  Ri  will  be  a  matrix  of  dimension  ni  x  no  where  "ni"  is  the  number  of 

inputs  and  "no"  is  the  number  of  outputs. 

To  obtain  the  pole-zero  form,  the  m  terms  in  the  pole-residue  form  above 
are  multiplied  out  to  form  a  ratio  of  polynomials  in  s  which  has 
denominator  of  order  2m  and  numerator  of  order  2(m-l).  The  numerator 
polynomial  can  then  be  factored  to  yield  the  zeros.  The  poles  are  already 
known  as  the  natural  frequencies  and  thus  the  denominator  does  not  have 
to  be  factored. 

A  state  space  realization  of  the  modal  equations  can  be  easily  assembled 
using  the  following  procedure.  The  matrix  equation  for  the  system  can  be 
rewritten  as  n  scalar  equations 

2  T 

rj .  +  to .  r|.=P  Gu/m.,  i=i . n. 

The  product  of  PjT  and  G  is  found  as  explained  above.  Defining  PjTG  /mj 
=Gmi  we  can  then  transform  this  second  order  system  to  the  first  order 
system, 
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n  =  Fti  4-  Gmu  , 

with  the  even  rows  of  Gm  being  composed  of  Gm.  and  the  odd  rows  of  Gm 
being  rows  of  zeros.  F  is  a  block  diagonal  matrix  with  the  i'th  sub  block 

ro  1 1 


The  output  distribution  matrix  for  the  state  space  realization  is 

y=nmn 

Hm=HP 

The  calculation  of  Hm  is  similar  to  the  calculation  of  Gm  but  has,  in 
general,  even  columns  composed  of  all  zeros  since  the  modal 
displacements,  and  not  velocities,  are  usually  output.  TRANSGEN  outputs 
the  second  order  state  space  form  since  this  requires  less  space,  so  the  user 
must  manually  perform  the  insertion  of  rows  or  columns  of  zeros  in  the  Gm 

and  Hm  matrices  and  the  expansion  of  the  Fj  sub-blocks.  This  is  easily  done 
using  common  matrix  manipulation  programs. 
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Examples 

Example  1.  Free-Free  Beam 

The  First  example  shows  how  TRANSGEN  would  be  used  to  Find  the 
diFferent  transfer  function  forms  for  a  continuous  uniform  beam  with  a 
force  actuator  perpendicular  to  the  beam  at  one  end  and  four  sensors  as 
shown  below  in  figure  Al.  The  beam  is  made  of  Aluminum,  has  a  length  of 
10  feet  and  square  cross  section  of  side  length  1  inch.  Sensors  1  and  3  are 
translation  sensors  while  2  and  4  are  rotation  sensors.  The  PAL  input  file 
in  Table  Al  shows  the  nodal  geometry  and  which  degrees  of  freedom  are 
made  active,  namely  z  translation  and  y  rotation. 

TITLE  NINE  POINT  BEAM 
NODAL  POINT  LOCATIONS  1 
10  00THROUGH9  12000 
10001 
1112001 

MATERIAL  PROPERTIES  1E7  0  2.6E-4  .3 
C  ALUMINUM  BEAM  OF  CROSS  SECTION  1  X  1  INCH 
BEAM  TYPE  211 

DO  CONNECT  1 2  THROUGH  8  9  STEP  1  1 

SPRING  IE-7 

CO  1  TO  10 

CO  9  TO  11 

ZERO 

TY  ALL 

TXALL 

RZ  ALL 

RXALL 

ALL  10,11 

ACTIVATE 

TZALL 

RYALL 

END  DEFINITION 

Table  Al.  PAL  Input  File 

The  reader  may  note  that  there  are  extra  springs  inserted  at  the  ends  of 
the  beam  perpendicular  to  it's  length.  These  very  small  springs  were  used 
to  force  PAL  to  obtain  rigid  body  mode  shapes  which  corresponded  to  pure  z- 
translation  and  pure  y-rotation  about  the  center  of  the  bar,  and  were  small 
enough  so  that  they  would  have  only  negligible  effects  on  the  flexible  mode 
shapes  and  frequencies.  Without  the  springs  PAL  returns  rigid  body  mode 
shapes  which  are  both  y-rotations  about  points  other  than  the  center  of 
mass  of  the  beam. 
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Tables  A3  and  A4  show  the  command  files  for  the  Normal  Modes  and 
Transient  Analysis  sections  of  PAL  dynamic  analysis  which  were  used  to 
obtain  the  mode  shape  and  modal  mass  files  seen  in  the.Tables  A6  and  A7. 


TITLE  FREEBEAM.MASS 
USE  FIRST  10  MODES 
JACOBI 

TIME  RANGE  1,1,1 

PRINT  MODAL  EQUATIONS 

SOLVE 

Table  A2,  PAL  Normal  Modes  Analysis  Command  File 


TITLE  FREEBEAM.DYN 
JACOBI 

RIGID-BODY  MODES  2 
SOLVE 

Table.  A3,  PAL  Transient  , Response  Analysis  Command  ..File 


Table  A4  below  shows  the  input  file  for  TRANSGEN  which  will  use  the 
first  eight  modes  (2  rigid  +  6  flexible)  to  calculate  the  desired  transfer 
functions  and  state  space  form  of  the  modal  equations.  Note  the 
correspondence  of  the  sensor  locations  specified  in  the  file  with  those 
pictured  in  figure  Al.  Also  note  that  even  though  there  are  only  nine  nodes 
on  the  beam  itself,  11  nodes  are  specified  in  line  three  since  there  are  two 
extra  nodes  (rigid  supports  for  the  small  springs).  This  file  will  call  for  all 
three  types  of  output  to  be  printed  in  the  output  file  which  the  user  is  asked 
to  identify  during  the  running  of  TRANSGEN. 
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freebeam.mode 

freebeam.mass 

11,8 

1.4 
1,3 

1.3 

1.5 

9.3 

9.5 

1,1,1 

Table  A4,  TRANSGEN  Input  File 


The  resulting  output  of  TRANSGEN  is  shown  in  Table  A8.  The  first 
section  of  the  output  file  contain  an  echo  of  the  relevant  input  information. 

Next  is  a  listing  of  the  zeros,  or  roots,  of  the  numerator  polynomials  for  all 
the  different  sensor  actuator  combinations.  We  thus  see  that  the  transfer 
function  from  actuator  1  to  sensor  1,  which  is  a  collocated  combination, 
has  the  following  form; 

y/s)  907,  Is2  (s2+60.62Hs2+196.  82)(s2  +411, 72)(s2+707.  52  )  (s2+1088 . 32  ) 
ui<s)  s2  s2  (s2  +  87.  94  2)  (s2  +  242.  52  )  (s2  +  475.  92 )  (s2  +  7882  )  (s2  +  1 1 832  ) . 

There  is  a  pole  zero  cancellation  at  the  origin  (the  two  zeros  cancel  two 
of  the  four  poles  at  the  origin) .  The  poles  and  zeros  then  alternate  with 
increasing  frequency  as  expected  for  a  collocated  sensor-actuator  pair.  The 
transfer  function  from  the  input  to  output  2  (rotation  sensor  at  left  end  of 
beam)  is  similar  but  with  different  zeros. 

The  transfer  functions  from  the  input  to  outputs  at  the  other  end  of  the 
beam  show  zeros  which  are  either  real  or  complex  instead  of  purely 
imaginary,  since  they  correspond  to  non-collocated  sensor  actuator  pairs. 

The  lower  zeros,  i.e.  those  at  8=  ±  122.4,  s=  ±307.3  are  quite  close  to  the  exact 
zeros  (exact  meaning  roots  of  the  exact  transfer  function  for  a  simple 
beam)  for  this  case  as  calculated  by  Wie  [A2],  The  higher  complex  zeros 
could  be  actual  zeros  or  they  could  be  the  results  of  numerical  truncation 
in  the  factoring  of  the  numerator  polynomial  performed  in  TRANSGEN. 

In  order  to  determine  if  they  are  actual  zeros  the  user  could  either  do 
another  TRANSGEN  run  with  more  modes  used  in  the  calculation  or  he 
could  take  the  state  space  description  output  from  TRANSGEN  and  use  a 
program  like  PC-  MATLAB  to  form  the  transfer  functions  independently. 

In  the  end,  it  will  most  likely  be  concluded  that  some  of  these  complex 
zeros  are  due  to  inaccuracies  of  the  mode  shapes  as  calculated  by  PAL,  and 
that  it  may  be  best  to  ignore  them  as  long  as  their  frequencies  appear  to  be 
much  higher  than  the  first  few  mode  frequencies. 

The  following  section  shows  the  natural  frequencies  of  the  beam  in 
units  of  radians/sec.  These  frequencies  are  really  the  imaginary  parts  of 
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the  roots  of  the  denominator  of  the  transfer  function,  also  known  as  poles. 
Each  pole  also  has  a  complex  conjugate  which  is  not  shown  here,  so  for 
this  example  there  will  be  four  poles  at  the  origin. 

The  next  section  of  output  data  is  the  modal  residues  section.  The 
residue  for  each  input-output  pair  is  shown  for  each  mode.  This  is  the  data 
that  is  used  to  compute  the  zeros  of  the  transfer  functions  as  explained  in 
the  theory  section  above.  The  residues  for  actuator  1  to  sensor  2  are 
combined  to  form  the  following  equation  which  is  then  multiplied  out  and 
factored  to  get  the  zeros. 

^ 2  _  0  1. 6  .  4.97  ,  8.41  11.83  ,  15.4 

ui  s2  +  88  s2  +  242.  5  s2  +  475.  9  s2  +  788 

Note  that  the  first  mode  has  a  residue  of  zero  since  the  rigid  translation 
mode  will  not  contribute  to  the  output  of  a  rotation  sensor  at  the  end  of  the 
beam.  The  other  residues  are  of  the  same  order  of  magnitude  here  only  by 
coincidence. 

The  last  section  of  output  data  is  the  state  space  form  of  the  modal 
equations.  The  G  matrix  is  8x1  since  there  is  one  input  and  eight  modes, 
while  the  H  matrix  is  4x8  since  there  are  four  outputs.  From  the  G  matrix 
the  user  can  note  that  all  8  modes  are  approximately  equally  controllable 
with  a  force  actuator  at  one  end  of  the  beam.  From  H  he  can  see  that  the 
first  mode  is  unobservable  with  the  rotation  sensors  (outputs  2  and  4)  as 
expected  since  it  is  a  rigid  body  translation.  The  G  and  H  matrices  can  be 
used  in  a  program  like  PC-MATLAB  or  similar  control  system  design 
software  after  insertion  of  rows  and  columns  of  zeros  as  described  in  the 
theory  section.  Natural  damping  inherent  to  all  structures  can  be  inserted 
in  the  F  matrix  if  the  analyst  already  knows  how  much  damping  each 
mode  contains. 

A  set  of  modal  state  space  equations  of  smaller  order  than  the  set 
which  is  output  can  be  formed  either  by  running  TRANSGEN  again 
requesting  a  smaller  number  of  modes  be  used,  or  by  simply  truncating 
the  larger  system  so  that  only  the  desired  number  of  modes  are  kept. 
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Table  A.1.  PAL  Normal  Modes  Analysis  Output 

FREEBEAM.DYN 

MCDE  SHAPE  NO.  1  AT  0.00000E-01  CPS  (  0.00000E-01  RAD/SEC) 


NODE 

X  TRANS 

Y  TRANS 

2  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0 . 0000E-C1 

0. OOOOE-Ol 

3 . 3331E-01 

0. OOOOE-Ol 

-3.7939E-07 

0 . CC00E-C1 

2 

0.0000E-01 

0. OOOOE-Ol 

3 . 3332E-01 

0. OOOOE-Ol 

-3.7939E-07 

0 . OOCOE-Qi 

3 

0 . 0000E-01 

0. OOOOE-Ol 

3 . 3332E-01 

0. OOOOE-Ol 

-3.7938E-07 

0 . OOOOE-Ol 

4 

0 . 0000E-01 

0. OOOOE-Ol 

3 . 3333E-01 

0. OOOOE-Ol 

-3 . 7936E-07 

o . occoe-o: 

5 

0 . 0000E-01 

0. OOOOE-Ol 

3 . 3333E-Q1 

0. OOOOE-Ol 

-3.7934E-07 

0. OOOOE-Ol 

6 

0 . 0000E-01 

0. OOOOE-Ol 

3 . 3334E-01 

0. OOOOE-Ol 

-3.7932E-07 

0. OOOOE-Ol 

7 

0 . 0000E-01 

0. OOOOE-Ol 

3 , 3334E-01 

0. OOOOE-Ol 

-3 . 7931E-07 

0. OOOOE-Ol 

8 

0 . OOOOE-Ol 

0. OOOOE-Ol 

3 . 3335E-01 

0. OOOOE-Ol 

-3.7930E-07 

0. OOOOE-Ol 

9 

0 . OOOOE-Ol 

0.  OOOOE-Ol 

3 . 3336E-01 

0. OOOOE-Ol 

-3.7929E-07 

0. OOOOE-Ol 

10 

0. 0000E-01 

0. OOOOE-Ol 

0.  OOOOE-Ol 

0.  OOOOE-Ol 

0.  OOOOE-Ol 

0. OOOOE-Ol 

11 

0. 0000E-01 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0.  OOOOE-Ol 

0. OOOOE-Ol 

MODE 

SHAPE  NO. 

2  AT  0 . 00000E-01  CPS  ( 

0 . 0Q000E-01  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0. 0000E-01 

0. OOOOE-Ol 

5 . 1624E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

2 

0. 0000E-01 

0. OOOOE-Ol 

3.8718E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

3 

0. 0000E-01 

0. OOOOE-Ol 

2 . 5812E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

4 

0. OOOOE-Ol 

0. OOOOE-Ol 

1 . 2907E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

5 

0. 0000E-01 

0. OOOOE-Ol 

1 . 1751E-05 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

6 

0. OOOOE-Ol 

0. OOOOE-Ol 

-1 . 2904E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

7 

0.  OOOOE-Ol 

0. OOOOE-Ol 

-2 . 5810E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

8 

0. OOOOE-Ol 

0. OOOOE-Ol 

-3 . 8716E-01 

0.  OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

9 

0. OOOOE-Ol 

0. OOOOE-Ol 

-5 . 1621E-01 

0. OOOOE-Ol 

8 . 6038E-03 

0. OOOOE-Ol 

10 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

11 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

MODE 

SHAPE  NO. 

3  AT  1 . 39966E+01  CPS  ( 

8 . 79434E+01  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0. OOOOE-Ol 

0. OOOOE-Ol 

5 . 5883E-01 

0. OOOOE-Ol 

2 . 1645E-02 

0. OOOOE-Ol 

2 

0. OOOOE-Ol 

0. OOOOE-Ol 

2 . 3668E-01 

0. OOOOE-Ol 

2 . 0994E-02 

0. OOOOE-Ol 

3 

0.  OOOOE-Ol 

0. OOOOE-Ol 

-5 . 5472E-02 

0. OOOOE-Ol 

1 . 7322E-02 

0. OOOOE-Ol 

4 

0. OOOOE-Ol 

0. OOOOE-Ol 

-2 . 6384E-01 

0. OOOOE-Ol 

9 . 8905E-03 

0. OOOOE-Ol 

5 

0. OOOOE-Ol 

0. OOOOE-Ol 

-3.3969E-01 

0. OOOOE-Ol 

-1 . 0394E-15 

0. OOOOE-Ol 

6 

0. OOOOE-Ol 

0. OOOOE-Ol 

-2 . 6384E-01 

0. OOOOE-Ol 

-9 . 8905E-03 

0. OOOOE-Ol 

7 

0. OOOOE-Ol 

0. OOOOE-Ol 

-5 . 5472E-02 

0. OOOOE-Ol 

-1.7322E-02 

0. OOOOE-Ol 

8 

0. OOOOE-Ol 

0. OOOOE-Ol 

2 . 3668E-01 

0. OOOOE-Ol 

-2 . 0994E-0? 

0. OOOOE-Ol 

9 

0.  OOOOE-Ol 

0. OOOOE-Ol 

5.5883E-01 

0. OOOOE-Ol 

-2 . 1645E-02 

0. OOOOE-Ol 

10 

0.  OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

11 

0.  OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 

0. OOOOE-Ol 
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mode 

SHAPE  NO. 

4  AT  3. 85917E+01  CPS  ( 

2.42479E+02  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0 . OQOOE-Ol 

0.0000E-01 

5.4647E-01 

0.0000E-01 

3 . 5812E-02 

0.0000E-01 

2 

0 . 0C00E-01 

0 . 0000E-01 

2. 6410E-02 

0.0000E-01 

3 . 1536E-02 

0 . 000CE-01 

3 

0 .  CC0CE-01 

C . 0000E-01 

-3.1972E-01 

0 . 0000E-01 

1 . 2231E-02 

0 . 00CCE-C1 

4 

0 . 00CCE-01 

0 . Q000E-01 

-3 . 0911E-01 

0 . 0000E-01 

-1 . 2990E-02 

0 . 0000E-C1 

5 

0 .  C000E-01 

0 . 0000E-01 

5.2976E-15 

0 . 0Q00E-01 

-2.4597E-02 

0 . 0000E-C1 

6 

0 . 0000E-01 

0.0000E-01 

3.0911E-01 

0 . 0000E-01 

-1 . 2990E-02 

0 . 0000E-C1 

7 

0 . 0000E-01 

0.0000E-01 

3.1972E-01 

0.0000E-01 

1 . 2231E-02 

O.OOOOE-Ol 

8 

0 . 0000E-01 

0 . 0000E-01 

-2. 6410E-02 

0 . 0000E-01 

3 . 1536E-02 

0 . 0C00E-C1 

3 

0 . 0000E-01 

0.0000E-01 

-5 . 4647E-01 

0 . 00G0E-01 

3 . 5812E-02 

O.OOOOE-Ol 

10 

0 . 0000E-01 

0.0000E-01 

0.0000E-01 

0.0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

11 

0 . 0000E-01 

0.0000E-01 

0.0000E-01 

0 . 0000E-01 

0.0000E-01 

O.OOOOE-Ol 

MODE 

SHAPE  NO. 

5  AT  7 . 57386E+01  CPS  ( 

4 . 75880E+02  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0 . 0000E-01 

0 . 0000E-01 

-5 . 3454E-01 

0.0000E-01 

-4 . 9055E-02 

O.OOOOE-Ol 

2 

0 . 0000E-01 

0 . 0000E-01 

1 . 4295E-01 

0 . 0000E-01 

-3 . 4982E-02 

O.OOOOE-Ol 

3 

0 . 0000E-01 

0 . 0000E-01 

3 . 3218E-01 

0 . 0000E-01 

1 . 1731E-02 

O.OOOOE-Ol 

4 

0 . 0000E-01 

0 . 0000E-01 

-7 . 8635E-02 

0 . 0000E-01 

3 . 3610E-02 

O.OOOOE-Ol 

5 

0 . 0000E-01 

0 . 0000E-01 

-3 . 8056E-01 

0 . 0000E-01 

0.0000E-01 

O.OOOOE-Ol 

6 

0 . 0000E-01 

0 . 0000E-01 

-7 , 8635E-02 

0.0000E-01 

-3 . 3610E-02 

O.OOOOE-Ol 

7 

0 . 0000E-01 

0 . 0000E-01 

3. 3218E-01 

O.OOOOE-Ol 

-1 . 1731E-02 

O.OOOOE-Ol 

8 

O.OGOOE-Ol 

0 . 0000E-01 

1 . 4295E-01 

0 . 0000E-01 

3.4982E-02 

O.OOOOE-Ol 

9 

0 . 0000E-01 

0 . 0000E-01 

-5.3454E-01 

0.0000E-01 

4 . 9055E-02 

O.OOOOE-Ol 

10 

0 . 0000E-01 

0 . 0000E-01 

0 . OOOQE-Ol 

0.0000E-01 

0.0000E-01 

O.OOOOE-Ol 

11 

0 . 0000E-01 

0 . 0000E-01 

O.OOOOE-Ol 

0.0000E-01 

0.0000E-01 

O.OOOOE-Ol 

MODE 

SHAPE  NO. 

6  AT  1 . 25510E+02  CPS  ( 

7.88600E+02  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0.0000E-01 

0.0000E-01 

-5 . 2213E-01 

0.0000E-01 

-6 . 1748E-02 

O.OOOOE-Ol 

2 

0 . 0000E-01 

0 . 0000E-01 

2 . 6383E-01 

0.0000E-01 

-2 . 9315E-02 

O.OOOOE-Ol 

3 

0 . 0000E-01 

0.0000E-01 

1 . 3305E-01 

0.0000E-01 

3 . 9374E-02 

O.OOOOE-Ol 

4 

0 . 0000E-01 

0.0000E-01 

-3.6447E-01 

0.0000E-01 

8.2770E-03 

O.OOOOE-Ol 

5 

0 . 0000E-01 

0 . 0000E-01 

0 . 0000E-01 

0.0000E-01 

-4 . 3599E-02 

O.OOOOE-Ol 

6 

0. 0000E-01 

0.0000E-01 

3 . 6447E-01 

0.0000E-01 

8 . 2770E-03 

O.OOOOE-Ol 

7 

0 . 0000E-01 

0 . 0000E-01 

-1 . 3305E-01 

0.0000E-01 

3 . 9374E-02 

O.OOOOE-Ol 

8 

0 . 0000E-01 

0 . 0000E-01 

-2 . 6383E-01 

0.0000E-01 

-2 . 9315E-02 

O.OOOOE-Ol 

9 

0 . 0000E-01 

0 . 0000E-01 

5.2213E-01 

0.0000E-01 

-6 . 1748E-02 

O.OOOOE-Ol 

10 

0 . 0000E-01 

0.0000E-01 

0 . 0000E-01 

0.0000E-01 

0.0000E-01 

O.OOOOE-Ol 

11 

0 . 0000E-01 

0 . 0000E-01 

0.0000E-01 

0.0000E-01 

0.0000E-01 

O.OOOOE-Ol 

MODE 

SHAPE  NO. 

7  AT  1 . 88251E+02  CPS  ( 

1 . 18282E+03  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 
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1 

0 . 0000E-01 

0 . 0000E-01 

-5.0768E-01 

O.OOOOE-Ol 

-7.3752E-02 

0. 2C0CE-C1 

2 

0 . 0000E-01 

0 . 0000E-01 

3 . 2585E-01 

O.OOOOE-Ol 

-1. 3803E-02 

0 . C000E-C1 

3 

0.0000E-01 

0 . 0000E-01 

-1 . 4433E-01 

O.OOOOE-Ol 

4.7186E-02 

0 . COOOE-Q1 

4 

0 . 0000E-01 

0 . OOOOE-Ol 

-2 . 0030E-01 

O.OOOOE-Ol 

-4 . 3204E-02 

O.OOOOE-Ol 

5 

0 . 0000E-01 

0.0000E-01 

3. 6168E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

6 

0 . 0000E-01 

0 . 0000E-01 

-2 . 0030E-01 

O.OOOOE-Ol 

4. 3204E-02 

O.OOOOE-Ol 

7 

0 . 0000E-01 

0 . 0000E-01 

-1.4433E-01 

O.OOOOE-Ol 

-4 . 7186E-02 

O.OOOOE-Ol 

8 

0 . 0000E-01 

0.0000E-01 

3 . 2585E-01 

O.OOOOE-Ol 

1 . 38Q3E-02 

O.OOOOE-Ol 

9 

0 . 0000E-01 

0 . 0000E-01 

-5 . 0768E-01 

O.OOOOE-Ol 

7.3752E-02 

0 . 000CE-01 

10 

0 . 0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

0 . 000CE-01 

11 

0 . 0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

*  — > 

L>  Z. 

SHAPE  NO. 

8  AT  2.641 

02E+02  CPS  ( 

1 . 65940E+03  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  RCT 

1 

0 . 0000E-01 

0.0000E-01 

4 . 8558E-01 

O.OOOOE-Ol 

8.4162E-02 

O.OOOOE-Ol 

2 

0 .  OOoOE-Ol 

O.OOOOE-Ol 

-3 . 2352E-01 

O.OOOOE-Ol 

-9 . 7318E-03 

O.OOOOE-Ol 

3 

0 . 0000E-01 

0.0000E-01 

3 . 2950E-01 

O.OOOOE-Ol 

-2 . 0942E-02 

O.OOOOE-Ol 

4 

0 . 0000E-01 

0.0000E-01 

-1.9839E-01 

O.OOOOE-Ol 

4 . 8086E-02 

O.OOOOE-Ol 

5 

0 . 0000E-01 

0.0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

-5 . 8190E-02 

O.OOOOE-Ol 

6 

0 . 0000E-01 

0.0000E-01 

1.9839E-01 

O.OOOOE-Ol 

4 . 8086E-02 

O.OOOOE-Ol 

7 

0 . 0000E-01 

O.OOOOE-Ol 

-3 . 2950E-01 

O.OOOOE-Ol 

-2 . 0942E-02 

O.OOOOE-Ol 

8 

0 . 0000E-01 

0.0000E-01 

3 . 2352E-01 

O.OOOOE-Ol 

-9.7318E-03 

O.OOOOE-Ol 

9 

0 . 0000E-01 

0.0000E-01 

-4 . 8558E-01 

O.OOOOE-Ol 

8 . 4162E-02 

O.OOOOE-Ol 

10 

0 . 0000E-01 

0.0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

11 

0 . 0000E-01 

0.0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

MODE 

SHAPE  NO. 

9  AT  3 . 49371E+02  CPS  ( 

2 . 19516E+03  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0 . 0000E-01 

O.OOOOE-Ol 

4 . 1905E-01 

O.OOOOE-Ol 

8.4679E-02 

O.OOOOE-Ol 

2 

0 . OOOOE-Ol 

O.OOOOE-Ol 

-2. 4673E-01 

O.OOOOE-Ol 

-3 . 2592E-02 

O.OOOOE-Ol 

3 

0 . 0000E-01 

0.0000E-01 

3 . 0171E-01 

O.OOOOE-Ol 

2 . 4485E-02 

O.OOOOE-Ol 

4 

0.0000E-01 

0.0000E-01 

-3.2716E-01 

O.OOOOE-Ol 

-1 . 2500E-02 

O.OOOOE-Ol 

5 

0.0000E-01 

O.OOOOE-Ol 

3 . 3608E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

6 

0 . 0000E-01 

O.OOOOE-Ol 

-3 . 2716E-01 

O.OOOOE-Ol 

1 . 2500E-02 

O.OOOOE-Ol 

7 

0 . 0000E-01 

O.OOOOE-Ol 

3 . 0171E-01 

O.OOOOE-Ol 

-2 . 4485E-02 

O.OOOOE-Ol 

8 

0 . 0000E-01 

O.OOOOE-Ol 

-2.4673E-01 

O.OOOOE-Ol 

3 . 2592E-02 

O.OOOOE-Ol 

9 

0 . 0000E-01 

O.OOOOE-Ol 

4 . 1905E-01 

O.OOOOE-Ol 

-8 . 4679E-02 

O.OOOOE-Ol 

10 

0 . 0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

11 

0 . OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

MODE 

SHAPE  NO. 

10  AT  4 . 86785E+02  CPS  ( 

3 . 05856E+03  RAD/SEC) 

NODE 

X  TRANS 

Y  TRANS 

Z  TRANS 

X  ROT 

Y  ROT 

Z  ROT 

1 

0 . 0000E-01 

O.OOOOE-Ol 

-5 . 9611E-01 

O.OOOOE-Ol 

-1 . 4691E-01 

O.OOOOE-Ol 

2 

0 . 0000E-01 

O.OOOOE-Ol 

2 . 1177E-01 

O.OOOOE-Ol 

9 . 7273E-02 

O.OOOOE-Ol 

3 

0 . 0000E-01 

O.OOOOE-Ol 

-1 . 5635E-01 

O.OOOOE-Ol 

-1.1626E-01 

O.OOOOE-Ol 

4 

0 . 0000E-01 

O.OOOOE-Ol 

8 . 0076E-02 

O.OOOOE-Ol 

1 . 2634E-01 

O.OOOOE-Ol 

5 

0.0000E-01 

O.OOOOE-Ol 

O.OOOOE-Ol 

O.OOOOE-Ol 

-1 . 2980E-01 

O.OOOOE-Ol 
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6  0 . OOOOE-Ol 

7  0. OOOOE-Ol 

8  0. OOOOE-Ol 

9  0. 0000E-01 

10  0. OOOOE-Ol 

11  0. OOOOE-Ol 


0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 


-8.0076E-02 
1 . 5635E-01 
-2 . 1177E-01 
5 . 9611E-01 
0. OOOOE-Ol 
0. OOOOE-Ol 


0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 


1 . 2634E-01 
-1.1626E-01 
9.7273E-02 
-1.4691E-01 
0. OOOOE-Ol 
0. OOOOE-Ol 


0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 
0. OOOOE-Ol 


Table  A.2.  PAL  Transient  Analysis  Output 

FREEBEAM . MASS 


MODAL  VALUES  FOR  M  MATRIX 


3 . 467E-03 

1. 958E-03 

2 . 772E-03 

1 . 788E-03 

2 . 437E-03 

1 . 569E-03 

2 . 327E-03 

3. 082E-03 

2 . 217E-03 

2.096E-03 

MODAL  VALUES 

FOR  C  MATRIX 

0 . 000E-01 

0 . 000E-01 

O.OOOE-Ol 

0 .  OOOE-Ol 

0. OOOE-Ol 

0. 000E-01 

0. 000E-01 

0. 000E-01 

0. 000E-01 

0 . C00E-01 

MODAL  VALUES 

FOR  K  MATRIX 

2 . 222E-08 
2.739E+03 

5 . 330E-08 
4.923E+03 

1 . 884E+01 

7 . 560E+03 

1 . 368E+02 

2 . 883E+04 

5 . 021E+02 

1 . 304E+03 

Table  A.3.  TRANSGEN  Output  File 


The  following  is  data  for  the  PAL  input  files 
freebeam.mode  freebeam.mass 

There  are  1  actuators  and  4  sensors. 
Actuators  location, orientation  are 
1,3 

Sensors  location, orientation  are 
1,3  1,  5  9,  3  9,5 


Zeros  of  element  1  1  of  the  transfer  function  numerator  matrix 


numerator  multiplying  factor  is  907.1 


0 . 000000E+00+ 
-0 . 286830E-36+ 
0 . 405669E-15+ 
0 . 405669E-1O+ 


0.000000E+00* j 
0 . 000000E+00*  j 
0 . 606458E+02* j 
-0 . 606458E+02*  j 


possible  pole  zero  cancellation 
possible  pole  zero  cancellation 
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-0 . 748325E-10+ 
-0.748325E-10+ 
C.423430E-07+ 
0. 423430E-07+ 
-0 . 156208E-05+ 
-0.156208E-05+ 
0. 991068E-05+ 
0 . 991068E-05+ 
-0.133909E-04+ 
-0.133909E-04+ 


0 . 196769E+03*  j 
-0 . 196769E+03*  j 
0.411693E+03* j 
-0 . 411693E+03*  j 
0 . 707490E+03* j 
-0 . 707490E+03*  j 
0 . 108826E+04*  j 
-0 . 108826E+04*  j 
0 . 156100E+04* j 
-0 . 156100E+04*  j 


Zeros  of  element  2  1  of  the  transfer  function  numerator  matrix 


numerator  multiplying  factor  is  84.16 


0 . 000000E+00+ 
-C .  212906E-35+ 
0.148755E-16+ 
0.148755E-16+ 
-0 . 225091E-10+ 
-0 . 225091E-10+ 
0.239905E-07+ 
0 . 239905E-07+ 
-0.115668E-05+ 
-0 . 115668E-05+ 
0 . 838468E-05+ 
0 . 838468E-05+ 
-0 . 122520E-04+ 
-0 . 122520E-04+ 


0. 000000E+00* j 
0 . 000000E+00*  j 
0 . 391307E+02* j 
-0 . 391307E+02* j 
0 . 157639E+03* j 
-0 . 157639E+03* j 
0. 357889E+03* j 
-0 . 357889E+03*  j 
0.643106E+03* j 
-0 . 643106E+03* j 
0.101919E+04* j 
-0 . 101919E+04*  j 
0 . 149776E+04*  j 
-0.149776E+04* j 


possible 

possible 


pole  zero  cancellation 
pole  zero  cancellation 


Zeros  of  element  3  1  of  the  transfer  function  numerator  matrix 

numerator  multiplying  factor  is  -65.70 


0 . 000000E+00+ 
-0 . 415497E-37+ 
0.122458E+03+ 
-0.122458E+03+ 
0.307296E+03+ 
-0.307296E+03+ 
0.425389E+03+ 
0 . 425389E+03+ 
-0 . 425389E+03+ 
-0.425389E+03+ 
0.497451E+03+ 
0. 497451E+03+ 
-0 . 497451E+03+ 
-0.497451E+03+ 


0.000000E+00* j 
0 . OOOOOOE+OO*  j 
0.000000E+00* j 
0 . 000000E+00*  j 
O.OOOOOOE+OO* j 
0. 0Q0000E+00* j 
0 . 359202E+03* j 
-0 . 359202E+03* j 
0 . 359202E+03* j 
-0 . 359202E+03* j 
0.104043E+04* j 
-0.104043E+04* j 
0 . 104043E+04*  j 
-0 . 104043E+04*  j 


possible  pole 
possible  pole 
possible  pole 
possible  pole 
possible  pole 
possible  pole 


zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 


Zeros  of  element  4  1  of  the  transfer  function  numerator  matrix 


numerator  multiplying  factor  is  12.34 


0.000000E+00+ 
-0.312078E-36+ 
0.787205E+02+ 
-0.787205E+02+ 
0 . 238CC5E<-03+ 
-0 . 238805E+03+ 


O.OOOOOOE+OO* j 
O.OOOOOOE+OO* j 
0. 000000E+00* j 
0.000000E+00* j 
0 . 000000E+00*  j 
0. 0Q0000E+00* j 


possible  pole 
possible  pole 
possible  pole 
possible  pole 
possible  pole 
possible  pole 


zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 
zero  cancellation 
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0.361175E+03+ 
0. 361175E+03+ 
-0.361175E+03+ 
-0 . 361175E+03+ 
0 . 454174E+03+ 
0 . 454174E+03+ 
-0. 454I74E+03+ 
-C.454174E+03+ 


3 . 338212E+03*  j 
-0 . 338212E+03*  j 
0 . 338212E+03*  j 
-0 . 338212E+03*  j 
0 . 101174E+04*  j 
-0.101174E+04* j 
0 . 101174E+Q4*  j 
-0.101174E+04* j 


THE  DENOMINTOR  ROOTS  ARE  (NATURAL  FREQUENCIES  rad/ sec) 

.0000 

.0000 

87.94 

242.5 
475.9 

788.6 


1659 . 


residue  matrix  for  mode  1  is... 
0. 3204E+02 
-0. 3647E-04 
0 . 3205E+02 
-0. 3646E-04 

residue  matrix  for  mode  2  is... 
0 . 9614E+02 
0 . 1602E+01 
-0 . 9614E+02 
0 . 1602E+01 

residue  matrix  for  mode  3  is... 
0 . 1281E+03 
0. 4963E+01 
0 . 1281E+03 
-0.4963E+01 


residue  matrix  for  mode  4  is... 
0 . 1283E+03 
0 . 8410E+01 
-0 . 1283E+03 
0 . 8410E+01 


residue  matrix  for  mode  5  is... 
0 . 1289E+03 
0.1183E+02 
0.1289E+03 
-0.1193E+02 


residue  matrix  for  mode  6  is... 
0 . 1301E+03 
0.1538E+02 
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-0 . 1301E+03 
0. 1538E+02 


residue  matrix  for  mode  7  is... 
0 . 1316E+03 
0 . 1912E+02 
0. 1316E+03 
-0 . 1912E+02 

residue  matrix  for  mode  8  is... 
0. 1319E+03 
0 . 2286E+02 
-0.1319E+03 
0 . 2286E+02 


G  matrix  is 
0 . 9614E+02 
0.1862E+03 
0 . 2293E+03 
0.2348E+03 
-0.2411E+03 
-0 . 2491E+03 
-0 . 2593E+03 
0.2716E+03 


H  matrix  TRANSPOSED  is 


0 . 3333E+00 
0 . 5162E+00 
0 . 5588E+00 
G . 5465E+00 
-0 . 5345E+00 
-0. 5221E+00 
-0.5077E+00 
0 . 4856E+00 


-0.3794E-06 
0 . 8604E-02 
0 . 2164E-01 
0 . 3581E-01 
-0.4905E-01 
-0 . 6175E-01 
-0.7375E-01 
0.8416E-01 


0 . 3334E+00 
-0 . 5162E+00 
0. 5588E+00 
-0 . 5465E+00 
-0 . 5345E+00 
0 . 5221E+00 
-0.5077E+00 
-0. 4856E+00 


-0 . 3793E-06 
0 . 8604E-02 
-0 . 2164E-01 
0 . 3581E-01 
0 . 4905E-01 
-0 . 6175E-01 
0.7375E-01 
0 . 8416E-01 
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c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 


TRANSGEN 


■  U  Uul  —  T  ,  « 


This  is  a  program  to  read  output  files  from  MSC  PAL  which 
contains  eigenvectors  for  each  mode  (nodal  displacements) . 

The  information  is  reordered  and  transfer  functions  and  residue 
are  found  for  an  input  sensor-actuator  set. 

Other  information  needed  is  a  modal  mass  and  stiffness  file  r.am 
and  the  locations  and  orientations  of  the  sensors  and 
actuators  . 

This  version  can  handle  a  max  of  30  modes  for  a  model  with 
any  number  of  node  points. 

This  program  is  freeware  and  is  not  meant  to  be  sold  or  exchang 
for  any  value 


iversity 


Written  by  Fred  Wiesinger,  Star.f or 


INTEGER  MX, resnum, FLG, resout, transout, ghout, actnum, sennum 
PARAMETER (MX=30) 
parameter (ma=6) 

DOUBLE  PRECISION  D (0 :MX) , R (MX) , IMG (MX) 

real*8  freq(MX) , hold(100, 6) , num (ma,ma,MX) , xmass (MX) , xstif (MX) 
BEAL*8  C (MX) , f req2 (MX) , specres (ma, ma, MX) , ASUM, reshold (MX) 
integer  actnods (ma) , aetdofs (ma) , sennods (ma) , sendofs (ma) 
real*8  g (MX,ma) , h (ma,MX) 

character*20  fname, namest, mkname, namepal, atest, alpha, outf ile, 
1  echofile 
integer  CK 


write (*, 1313) 

1313  FORMAT (IX, 'THE  REQUIRED  INPUTS  ARE: ',/, 'NAME  OF  PAL  MODAL  OUTPUT 

1  FILE' ,/, 'NAME  OF  PAL  MODAL  MASS  AND  STIFFNESS  FILE  (FROM  TRANS.  AN 

2  .)',/, 'NUMBER  OF  NODES  IN  THE  MODEL',/, 'NUMBER  OF  MODES  TO  BE  USED 

3  ',/,'  ACTUATOR  AND  SENSOR  NODAL  PLACEMENT  AND  ORIENTATION',/, 

4  'ENTER  A  "1"  IF  YOU  WOULD  LIKE  TO  SPECIFY  THESE  BY  HAND  ',/, 

5  'ENTER  A  "2"  IF  THEY  SHOULD  BE  READ  FROM  A  FILE') 

READ  (*,  *)  PARM 

IF ( PARM= 1 )  THEN 

WRITE (*,*) 'ENTER  PAL  MODES  FILE  NAME  IN  QUOTES’ 

READ (*,  *)  namepal 

write  (*,*)  'ENTER  PAL  MODAL  MASS  FILE  NAME  IN  QUOTES' 

READ  ( * , * )  mkname 

1  write (*,*)  'ENTER  NUMBER  OF  NODES  IN  STRUCTURE  AND  NUMBER 
10F  MODES  TO  BE  USED' 

READ ( * , * )  nodnum,modenum 

WRITE (*,*)  'ENTER  NUMBER  OF  ACTUATORS,  NUMBER  OF  SENSORS  ' 

READ (*,  * )  actnum, sennum 
do (i=l, actnum) 

write(*,*)  'ENTER  ACTUATOR  #  ' , i, ' -  NODE  NUMBER  AND 
1  D.O.F. (l->  6) ' 

READ (*, *)  actnods (i) , aetdofs (i) 
repeat 

do (i=l, sennum) 
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write  (*,  *)  'ENTER  SENSOR  X  NODE  NUMBER  AND  D.C.F.  (l-> 

READ  ( * ,  * )  sennods (i) , sendofs (i) 

repeat 

write  (*,*)  'Enter  l="yes", 0="no"  for  the  following  output 
1  options' 

write(*,*)  'Residues  output, zeros  output, g  &  h  matrices  cutpu 
read!*, *)  resout, transout, ghout 


elseif (PARM=2 ) THEN 

WRITE (*,*)  'ENTER  INPUT  FILE  NAME' 
read(*,*)  fname 
open (unit=12, file=fname) 
c  read  name  of  PAL  output  file 
read(12, 1000)  namepal 
1000  format (A20) 

read (12, 1009)  mkname 

1009  format (a20) 

c  read  number  of  Nodes,  and  number  of  Modes  to  be  used.**** 
read (12,*)  nodnum,modenum 
c  read  number  of  actuators  and  sensors  **** 
read (12,*)  actnum, sennum 
do  (i=l, actnum) 

c  read  actuator  nodes  and  d.o.f.'s. 

READ (12,*)  actnods (i) , actdofs (i) 
repeat 

c  read  sensor  nodes  and  d.o.f.'s. 
do (i=l, sennum) 

READ (12,*)  sennods (i) , sendofs (i) 
repeat 

c  Enter  l="yes", 0="no"  for  the  following  output  options 
c  Residues  output, zeros  output, g  &  h  matrices  outDut 
READ (12,*)  resout, transout, ghout 

close (unit=12) 

else 

GO  TO  999  !  go  to  end 

endif 

c  Input  modal  mass  and  stiffness  from  PAL  output  file  for  each  mode 
c  Modal  mass  and  stiffness  are  generated  by  PAL's  transient 

c  analysis  when  the  "print  modal  equations"  command  is  given, 

write  (*, *)  mkname,'  IS  MASS  MATRIX  FILE' 
open (unit=14, file=mkname) 
read (14, 1010) 

1010  format (//////) 

C  READ  MODAL  MASS  FOR  EACH  MODE 

if (modenum. It. 6)  go  to  41 
do  40  i=l, modenum/ 6 

40  read(14, 1011)  (xmass (i*6+ j-6) , j=l, 6) 

1011  format (lx, 6el2 . 4) 

41  l=modenum-modenum/6*6 
if (l.eq.0)  go  to  48 
lk=modenum-l 

read(14, 1011)  (xmass (i) , i=lk+l, modenum) 
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48  continue 

C  FIND  START  OF  STIFFNESS  MATRIX 

do  50  i=l, 200 
read(14, 1012)  alpha 
1012  format (a3) 

atest='  M'  ! searces  for 

in  medal 

if (atest.eq. alpha)  go  to  70  Ithird  one  is  medal 

stiffness 

GO  TO  50  ! section 

70  alphtes=l+alphtes 

if (alphtes . ge . 2) go  to  59 
50  continue 
c 

C  READ  MODAL  STIFFNESS  FOR  EACH  MODE 

59  read(14,*) 

if (modenum. It. 6)  go  to  61 
do  60  i=l, modenum/ 6 

60  read(14, 1011)  (xstif (i*6+ j-6) , j=l, 6) 

61  l=*nodenum-modenum/6*6 
if (1 .eq. 0)  go  to  68 
lk=modenum-l 

read(14, 1011)  (xstif (i) , i=lk+l, modenum) 

68  continue 

close (unit=14) 

c 


write  (*,*)  'enter  an  output  file  name  ' 
read(*,*)  outfile 

open (unit=15, f ile=outf ile, status=' new' ) 
write  (15, 1111)  namepal ,  mkr.ame 
write (15, 1112) actnum, sennum 

write (15, 1113) (actnods (i) , aetdof s (i) , i=l, actnum) 
write (15, 1114) (sennods (i) , sendofs (i) , i=l, sennum) 

1111  format (lx, ' The  following  is  data  for  the  PAL  input  files',/ 
1, 4X,a20,a20) 

1112  format (lx, 'There  are  ',11,'  actuators  and  ',11,'  sensors.') 

1113  format (2x, 'Actuators  location, orientation  are',/,3x, 

li3, ', ' , il, 2x, i3, ', ' , il, 2x, 13, ', • ,il,2x,i3, ' , ' ,il, 2x,i3, ' , ', 
2il,  2x, 13, ', ',il,2x,i3, ', ',il,2x,i3, ', ' ,il,2x) 

1114  format (2x, ' Sensors  location, orientation  are',/,3x, 

li3, ', ' , il, 2x, i3, ', ' , il, 2x, i3, ', ' , il, 2x, i3, ', ' , il, 2x, i3, ', ', 
2il,2x,i3, ', ',il,2x,i3, ', ',il,2x,i3, ', ' , il, 2x) 


c 


open (unit=13, f ile=namepal) 

do  900  n=l, sennum 
do  899  m=l, actnum 

actnod=actnods (m) 
actdof=actdofs (m) 


*■  he  M'*"n 
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sennod=sennods (n) 
sendof=sendofs (n) 


c  ”**  begin  the  data  read  from  PAL  output  file.***’** 
rewind (unit=13) 
read (13, 1007) 

1007  format (//) 

c  For  each  mode  read  the  frequency  and  the  values  of  the  modal 
displacements 

c  in  the  corresponding  degree  of  freedom  directions  and  use  the 
displacements 

c  to  calculate  the  residue  for  the  given  sensor  and  actuator  transfer 
c  function. 

do  42  i=l,modenum 
read(13, 1001)  freq(i) 

1001  format (/, 45x, ell . 5, ///) 

do  5  j=l,nodnum 

if ( j=actnod. or . j=sennod)  then 
read  (13, 1002)  (hold(  j,  k) ,  k=l,  6) 
else 

read (13, *) 
endif 

1002  format (lx, 6x, ell . 4, lx, ell . 4, lx, ell . 4, lx, ell . 4, lx, ell . 4 , lx, ell . 4 ) 

5  continue 

c  ***  calculate  the  residue  for  each  mode  ***** 

num (n,m, i)  =  (hold (actnod,  actdof ) ) * (hold(sennod, sendof ) ) 
c  *****  NOTE:  the  g  and  h  matrices  which  are  output  are  the  Gm  and  Hn 
c  ****  matrices  of  the  theory  section  of  the  TRANSGEN  user  manual, 

g (i,m) =hold (actnod, actdof) /xmass (i) 
h (n, i) =hold(sennod, sendof) 

42  continue 

if (resout=0 . and. transout=0)  go  to  899 

c 


C  FORM  SPECIFIC  RESIDUES  (residue  divided  by  modal  mass  for  each  mode) 

C  If  a  specific  residue  is  equal  to  zero  THAT  WHOLE  MODE  WILL  BE 

C  OMMITED  FROM  THE  TRANSFER  FUNCTION  CALCULATION. 

lowres=0 

do  (i=l,modenum) 

specres (n,m, i) =num (n,m, i) /xmass (i) 
if (specres (n,m, i) =0 . )  lowres=lowres+l 
repeat 

if (lowres=modenum) then 

write (15,*)  ’all  residues  are  0, sensor  ' ,n, '  input  is  not 
lassociated  with  actuator  ',m, ’output’ 
go  to  899 
endif 

if (transout=0)go  to  899 


! define  freq2  and  reshrl; 


vector 


do  (i=l,modenum) 

reshold(i) =specres (n,m, i) 
freq2(i)=  -freq(i)**2 
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repeat 


resnum=modenum-lowres 
if(resnum=l)  then 

write  (*,*)  'only  one  residue  is  nonzero  so  t.f.  is  trivial' 
go  to  899 
endif 
c 


c  compute  products  of  partial  fractions  (i.e.  multiply  out  the 

partial 

c  fractions. 

call  pf (modenum, resnum, reshold, freq2, C, ASUM, FLG) 

c 


C  NUMERATOR  COEFICIENTS  OUTPUT  FROM  pf  ARE  ORDERED  AS  BELOW. . . 

C  ASUM*  { SAN+  C  (N)  *SAN-1+  C (N-l) *SAN-2+  ...  +C(2)*S+  C(l)} 

C  WHERE  N=modenum-l, ASUM=numerator  multiplying  factor 

C 

C  REORDER  COEFICIENTS  FOR  BAIRSTOW  ALGO. 

C  ODD  COEFICIENTS  ARE  ZERO  SINCE  SYSTEM  IS  UNDAMPED, N=2* (MODENUM-1 ) 
150  DO  (1=0, MX) 

D (I ) =0 . 

REPEAT 

IF ((FLG=1) .AND. (resnum. eq.l))  then 

write  (*,*)  'numerator  is  a  constant  =  ',ASUM 

go  to  899 

endif 

N=2* (resnum-1) 

D (0) =1 
DO  (1=1, N) 

CK=I/2 

IF (2*CK-I . LT . 0) THEN 

D (I) =.00001  ! BAIRSTOW  ALGO.  NEEDS  ALL 

COEFICIENTS  *■  0 
ELSE 

D  (I)  =C  ( (N-I+2)  /  2) 

ENDIF 

REPEAT 


c  call  subroutine  to  factor  transfer  function  numerator 
C  MASROOT  FINDS  ROOTS  OF  NUMERATOR  POLYNOMIAL 
call  MASROOT (N,D,R, IMG) 
c 


write  (15, *) 

write  (15,*)  'Zeros  of  element  ' ,n, '  ',m, 

1  '  of  the  transfer  function  numerator  matrix' 

write  (15,*)  'numerator  multiplying  factor  is  ',ASUM 
do (i=l, N) 

c  if (IMG (i) =-IMG (i+1) . and. i . ne .N)  go  to  1876 

f lagl=0 

do  (j=l, modenum) 
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if (abs ( (abs (IMG (i ) )-freq(j) ) / (freq( j) +le-7) ) .it.le-2)  flagl=l 
repeat 

if (flagl=l) then 

write  (15,  4041)  R(i),IMG(i) 

else 

write (15, 4040)  R(i),IMG(i) 
endif 

1376  repeat 

4040  format (4x, el4 . 6, ' +  ' , el4 . 6, ' * j ' ) 

4041  format (4x,el4 . 6,  '  +  ',el4.6,'*j  possible  pole 
1  zero  cancellation  ' ) 


899  continue 

900  continue 

close (unit=13) 

C  END  OF  MAIN  LOOPS 


write  (15, *) 

WRITE (15,*)  'THE  DENOMINTOR  ROOTS  ARE  (NATURAL  FREQUENC 
1IES  rad/ sec) ' 

DO  (I=l,modenum) 

WRITE (15,*)  freq(I) 

REPEAT 


if(resout=0)  go  to  1400 

c  print  residue  matrices  to  the  output  file  in  mode  order 
do  (i=l,modenum) 
write (15, 1555)  i 

1555  format (/, lx, ' residue  matrix  for  mode  ',i2,'  is...') 
do (n=l, sennum) 

write (15,1395)  (specres (n,m, i) ,m=l, actnum) 

repeat 

repeat 

write (15, 1398) 

1395  format (lx, 6(el2. 4) ) 

1398  format (///) 


1400  if(ghout=0)  go  to  999 

write  (15,*)  'G  matrix  is' 
do  (i=l,modenum) 

write (15, 1450) (g(i, j) ,  j=l, actnum) 

repeat 

write (15, *) 

write  (15,*)  'H  matrix  TRANSPOSED  is' 
do  (i=l,modenum) 

write  (15, 1450)  (h(j,i), j=l, sennum) 
repeat 

1450  format (2x, 6 (el2 . 4) ) 
c 
c 
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999  write (*,*)  'Program  complete.  Please  hit  return.’ 
pause 
stop 
end 

c  **********  SUBROUTINES  FOLLOW  ***************************** 


C  **********  p.p.  xo  TRANS  **************************  1/13/89  FW 
C  THIS  IS  A  SUBROUTINE  WHICH  WILL  FIND  THE  TRANSFER  FUNCTION  PCLE/ZERO 
C  FORM  FROM  A  POLE-RESIDUE  FORM. 

C  THE  INPUTS  ARE  N  (THE  NUMBER  OF  MODES ),Q  (THE  NUMBER  OF  NONZERO 
C  RESIDUES,  THE  N-VECTOR  A  (SPECIFIC  RESIDUES), 

C  AND  THE  N-VECTOR  SIG  (NEG  OF  MODAL  FREQUENCIES  SQUARED) . 

C  THIS  SUBROUTINE  CALLS  THE  SUBROUTINE  POLYMUL  TO  OBTAIN  NUMERATOR 
C  POLYNOMIAL  COEFICIENTS  FROM  NUMERATOR  ROOTS. 

C  INPUTS:  N=modenum,Q=resnum, A=specres, SIG=freq2 
C  OUTPUTS:  TOR, ASUM, FLG, (and  Q  could  be  updated) 
subroutine  pf (N, Q, A, SIG, TCR, ASUM,  FLG) 

INTEGER  Q,FLG 
PARAMETER (MX=30) 

REAL*8  A (MX) , SIG (MX) , RR (MX) , CR (MX) , TCR (MX)  ,ASUM 

C 

C  MULTIPLY  OUT  THE  PARTIAL  FRACTIONS,  THEN  ADD  POWERS  OF  S  COEFICIENTS 

C  AND  DIVIDE  EACH  COEFICIENT  BY  THE  COEFICIENT  OF  THE  HIGHEST  POWER  OF  S. 

ASUM=0 . 

DO  10  1=1, N 
L=0 

DO  20  J=1,N 
IF(J.EQ.I)  GO  TO  20 
IF (A ( J) =0 . )  GO  TO  20 
L=L+1 

RR  (L)  =SIG  ( J) 

20  CONTINUE 

C  INPUT  TO  POLYMUL  IS  VECTOR  RR,  AND  #  OF  FACTORS  TO  BE  MULTIPLIED  OUT. 

C  OUTPUT  IS  VECTOR  CR  WHICH  CONTAINS  POLYNOMIAL  COEFFICIENTS. 

CALL  POLYMUL (RR,Q-1,CR) 

DO  30  K=1,Q-1 
IF(I.EQ.l)  TCR(K) =0 

C  ADD  UP  ALL  THE  COEF . 'S  OF  LIKE  POWERS  OF  S  ... 

TCR (K) =TCR (K) +CR (K) *A ( I ) 

30  CONTINUE 

ASUM=A ( I ) +ASUM 
10  CONTINUE 

******  OUTPUT  ****** 

IF (ASUM=0 . )  THEN  ! IF  ASUM=0  SET  ASUM=NEXT  COEFICIENT  AND 

FLG=1  !  PROCEED 

ASUM=TCR(Q-1) 

Q=Q-1 
ENDIF 

DO  40  1=1, Q-l 
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TCR ( I ) =TCR ( I ) / ASUM  ! DIVIDE  EACH  COEFICIENT  BY  LEADING 

COEFICIENT 

40  CONTINUE 

close (unit=12) 

return 

END 

C 

C 

C  THIS  IS  A  SUBROUTINE  WHICH  WILL  RETURN  THE  COEFICIENTS  OF  A 
C  POLYNOMIAL  WHOSE  ROOTS  ARE  SENT  TO  IT. 

C  THE  INPUTS  ARE  THE  ROOTS  (I.E.  -5  FOR  (S+5)  ) , AND  THE  NUMBER  OF 
C  ROOTS,  N.  THE  OUTPUT  IS  AN  N  VECTOR  CR  CONTAINING  THE  POLYNOMIAL 
C  COEFICIENTS  AS  SHOWN  BELOW: 

C  P(S)=SAN  +CR(N) *SA (N-l)  + _ +C (2) *S  +C(1) 

0  ★★*★★★★★★★★★★★★★★★★★★★★★★★★★*★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★* 

SUBROUTINE  POLYMUL (RR,  N, CR) 

REAL* 8  CR(100) ,RR(100) ,DR(100) 

INTEGER  N 
IF(N.NE.l)  GO  TO  3 
CR(2)=1 
CR  (1)  =-RR  (1) 

2  RETURN 

3  DO  10  1=2, N 
CR(I) =0 

10  CONTINUE 
CR (2) =1 
CR (1 ) =-RR (1) 

DO  20  1=3, N+l 
DO  30  J=1,I-1 
DO  40  L=l, I 
DR (L) =CR (L) 

40  CONTINUE 

CR (I-J+l ) =DR (I- J) -DR (I-J+l) *RR (1-1 ) 

30  CONTINUE 
DR ( 1 ) =CR ( 1 ) 

CR(1) =-DR(l) *RR(I-1) 

20  CONTINUE 

DO  90  1=1, N 
WRITE (*, 1000)  CR(I) 

1000  FORMAT (IX, ’CR"S=  ' ,8F12.8) 

90  CONTINUE 
RETURN 
END 


C  ********  D (I) =  COEFFICIENT  OF  X**(N-I)  ************* 

C  Subroutine:  Finds  roots  of  polynomial  using  Bairstow  algorithm. 

C  Bairstow  works  fine  for  non-repeated  roots  and  non-zero  coeficients. 
C  N  is  order  of  polynomial 

C  D  must  be  of  dimension  N+l;  best  to  divide  through  by  D(0),  the  coef 
C  of  x**n  before  call  to  ROOTS 

C  Outputs:  IMG  is  imag  part  of  roots,  R  is  real  part  of  roots 
C  IERR=1  if  there  is  a  error  in  converging 
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SUBROUTINE  MASROOT (N, D, R, IMG) 

IMPLICIT  NONE 

INTEGER  I, J, J2, N, II, L, MX,  IERR 
PARAMETER (MX=30 ) 

DOUBLE  PRECISION  PI (MX) ,PR(MX) ,D(0:MX)  !  inputs 

DOUBLE  PRECISION  R (MX) , IMG (MX)  !  outputs 

DOUBLE  PRECISION  Q (0 :MX) , P (0 :MX) , B (0 :MX)  !  work  vars 
DOUBLE  PRECISION  E, R2, S2, R3, S3, R4, S4, R5, S5, D6, R6, S6 
DOUBLE  PRECISION  RA, RB, IA, IB 
IERR=0 
DO  (1=1, MX) 

R(I) =0 . 

IMG (I) =0 . 

REPEAT 
DO  (1=1,  N) 

PR(I) =0 . 

PI  (I)=0. 

REPEAT 

J=1 

J2=0 

E=. 0000001 
680  L=0 

IF  ( J.NE.N)  THEN 

R2=-PR(J) -PR( J+l) 

S2=PR ( J) *PR( J+l) +PI (J)  **2 
Il=N-2+J2 

IF  (I1.EQ.2)  GOTO  720 
IF  (I1.GT.2)  GOTO  730 
END  IF 
R  ( J)  — D(l) 

GOTO  910 
720  R2=D (1) 

S2=D (2) 

CALL  QUAD  (R2,  S2,  J,  RA,  RB,  IA,  IB) 

R  ( J)  =RA 
R ( J+l ) =RB 
IMG ( J) =IA 
IMG (J+l) =IB 
GOTO  910 
730  B (0) =1 

B (1) =D (1) -R2 
DO  (1=2,11) 

B (I) =D (I) -R2*B (I-l)-S2*B(I-2) 

REPEAT 
P(0)=0 
P(l)=-1 
DO  (1=2,11) 

P (I) =-B (1-1) -R2*P (1-1) -S2*P (1-2) 

REPEAT 
Q (0) =0 
Q(1)=0 
DO  (1=2,11) 

Q  (I) =-B (1-2) -R2+Q (1-1) -S2*Q (1-2) 

REPEAT 
R3=B ( II— 1 ) 

S3=B (II) +R2+B (II— 1 ) 
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R4=P  (II— 1 ) 

S4=P (II) +R2*P(I1-1) +B(I1-1) 

R5=Q(I1-1) 

S5=Q (II) +R2*Q (11-1) 

D6=R4*S5-R5*S4 
IF  (ABS (D6) .LT.E)  THEN 
GOTO  810 
ELSE 

R6=(S3*R5-S5*R3)/D6 

ENDIF 

S6=(R3*S4-R4*S3) /D6 

R2=R2+R6 

S2=S2+S6 

IF  ( (ABS (R6) . GT.E*ABS (R2) ) .OR.  (ABS ( S 6 ) . GT.E*ABS (S2) ) )  GOTO  820 
810  CALL  QUAD  (R2, S2,  J,  RA,  RB,  IA,  IB) 

R(J)=RA 
R (J+l ) =RB 
IMG  ( J)  =IA 
IMG ( J+l) =IB 
GOTO  850 

820  IF  (L.GE.N*50)  THEN 

WRITE (*,*)  '  Polynomial  solver  did  not  converge' 

IERR=1 
RETURN 
ENDIF 
L=L+1 
GOTO  730 
850  J2=J2+1 
J=J+2 

DO  (1=1,11) 

D ( I ) =B ( I ) 

REPEAT 
GOTO  680 
910  WRITE (*, * ) 

RETURN 

END 

*  Solves  quadratic  equation 

SUBROUTINE  QUAD (R2, S2, J, RA, RB, IA, IB) 

IMPLICIT  NONE 
INTEGER  J 

DOUBLE  PRECISION  DI, R2, S2, RA, RB, IA, IB 
R2=R2/2. 

DI=R2**2-S2 
IF  (DI.GE.0.)  THEN 
DI=SQRT(DI) 

RA-DI-R2 

RB=~DI-R2 

IA-0. 

IB=0 . 

ELSE 

DI=SQRT(-DI) 

RA=-R2 
RB=-R2 
IA=DI 
IB— DI 
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ENDIF 

RETURN 

END 
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